Kondo Phase Transitions of Magnetic Impurities in Carbon Nanotubes 
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We propose carbon nanotubes (CNTs) with magnetic impurities as a versatile platform to achieve 
unconventional Kondo physics, where the CNT bath is gapped by the spin-orbit interaction and 
surface curvature. While the strong-coupling phase is inaccessible for the special case of half-filled 
impurities in neutral armchair CNTs, the system in general can undergo quantum phase transi- 
tions to the Kondo ground state. The resultant position-specific phase diagrams are investigated 
upon variation of the CNT radius, chirality, and carrier doping, revealing several striking features, 
e.g., the existence of a maximal radius for nonarmchair CNTs to realize phase transitions, and an 
interference-induced suppression of the Kondo screening. We show that by tuning the Fermi energy 
via electrostatic gating, the quantum critical region can be experimentally accessed. 

PACS numbers: 73.22.-f, 73.20.Hb, 72.15.Qm, 64.70.Tg 
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Carbon nanotubes (CNTs) are formed by wrapping a 
graphene sheet into a cylinder of nanometer radius 
Their exceptional electronic structure has allowed the ex- 
ploration of various fascinating Kondo phenomena, in- 
cluding the singlet-triplet Kondo resonance [2J, the en- 
hanced shot noise due to the SU{4) Kondo effect 0], 
and the competitions with ferromagnetism [B[ as well as 
superconductivity [(J. These studies have utilized short 
CNTs to construct quantum dots behaving as artificial 
magnetic impurities. Long CNTs, on the other hand, can 
play the role of one-dimensional (ID) host for a real mag- 
netic impurity, which may be either a magnetic adatom 
on the top of a carbon atom (T site) or at the center of 
a hexagon (C site), or a substitutional dopant in a car- 
bon vacancy (S site). Indeed, the Kondo effect for cobalt 
clusters adsorbed on metallic CNTs has already been ob- 
served [7] . This has spurred several theoretical works 8 
Hfjj to address related issues. However, a generic Kondo 
model of a CNT-hosted magnetic impurity, pertaining 
to arbitrary positions at the atomic scale, has not yet 
been established. More importantly, while these theories 
have all considered the metallic-CNT host as ID normal 
metal to yield conventional Kond o physi cs recent 
experiments 12 1 and calculations l3-15| have demon- 
strated that metallic CNTs can not be normal metal, al- 
ways having rich low-energy band structures due to the 
spin-orbit interaction (SOI) and curvature effect. Inter- 
esting Kondo physics then arises when these ingredients 
in the CNT host are included. 

In this paper, after establishing a generic Hamiltonian 
for magnetic impurities in metallic CNTs, we show that 
depending on explicit impurity positions, the system can 
maps onto two kinds of hard-gap Kondo models whose 
host density of states (DOS) are identically gapped by the 
SOI and curvature effect, but scale distinctly outside the 
gap region due to the quantum interference between dif- 
ferent hybridization paths. We combine renormalization 
group (RG) arguments and slave boson (SB) techniques 
to demonstrate that the local-moment (LM) state per- 



sists for half-filled impurities in neutral armchair CNTs 
due to the particle-hole (p-h) symmetry. Away from this 
special case, quantum phase transitions (QPTs) exist in 
these gapped systems, separating the Kondo and LM 
ground states. The resultant phase diagrams are char- 
acterized by the CNT radius, chirality, carrier doping, 
and the impurity position. The interference is found to 
reduce the Kondo regime, making S/C configurations un- 
favorable for Kondo screening as compared with T sites. 
For sufficiently deep impurity levels, two quantum crit- 
ical regions are accessible by scanning tunneling probes 
and gating the CNT host, with signatures characterizing 
the nonarmchair from armchair chilarities. 

Our starting point is the Anderson Hamiltonian of a 
magnetic impurity in graphene, H—ri +'Hg+'H g, where 
W = Y^ a £ ddid a + ^Y^, a i\iai\ia models the impurity 



as usual 



rig is the graphene tight-binding Hamilto- 



nian reading H g = Y2(i,j), a to £ + H -C, here 
a CT (R a ,) [6o-(Rhj)] annihilates an 7r-band electron on sub- 
lattice A (B) at position R a i (Ry), and t is the nearest- 
neighbor hopping energy. The hybridization term H. g = 
£ CT gl d CT +H.c. with gl=Y, {j) [K;4(Rai)+V^(R«)] , 
where (J) stands for the A and/or B sublattice sites 
nearest to the impurity, and V x j (x=a,b) represent the 
corresponding hybridization amplitudes. In particular, 
gj T =Vaia' (T (0) for a T-site adatom. C-site impurities can 
hybridize with six surrounding carbon atoms, yielding 
<?£=Ei=i [Ktf4(Xj) + Vbjbt (-Xj-)] with X, the lattice 
nearest-neighbor vectors, while gt=Ylj=i Hj&o-ptj) for 
S'-site impurities on sublattice A. In momentum space, 
the fermionic basis c kaa =^{a.a\ icr J r-^-^b] icr ) diagonalizes 

the graphene Hamiltonian as Ug=J2k, a ,a £ a( k -)cl cra c\ i aa, 
where fa = J2j=i e lk ' X 3, e Q (k) = ai|</> k |, and a = ±1. 
Close to the Dirac points K, the dispersion is linear, 
i.e., e a (K + k) ~ ahvp\n\ for \n\ <C |K|, with vp the 
Fermi velocity. In this basis, the hybridization becomes 
Kog=T,u,*, a V a (k)cl aa d a +B.c, where y Q (k)=(a$ ak + 

and N 
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the number of sublattice sites. 

We now roll up the graphene sheet along the chiral vec- 
tor C/j=niai + ri2a.2 to create a (rti, 712) CNT [![, where 
711,712 € Z and ai,a2 are the primitive lattice vectors. 
While k's component parallel to the tube axis, p=K\i , re- 
mains continuous for CNTs of long length L, the periodic 
boundary condition, (K+k) •C/j=27rm, m € Z, quantizes 
k's perpendicular component, q T = kj_ = (m + tv/3)/R. 
Here i? = #r yTif + 7i| + 7ii7i2 is the tube radius with 
a = |ai| ~ 2.46A the lattice constant, the valley in- 
dex t = ±1 denotes the two inequivalent K + ,K _ Dirac 
points, and i/=mod(ni — 712, 3) characterizes the metal- 
lic (u = 0) or semiconducting (i/ = ±l) CNTs. Restrict- 
ing the graphene quantities, e Q (k), V Q (k), Cko-a, only to 
these allowed Bloch states near the Dirac points yields 
corresponding quantities for the CNT: the 7r-band spec- 

1 the hybridization V pTa — 



trum e„ TQ = aHvFy/p 2 



It 



V a (K T + k)\ k= ,., and the operator c paTa - 

The surface curvature of CNTs induces the it band 
hybridizing with other high-energy bands (CIH), and en- 
hances the effect of intrinsic SOI, V so , of carbon atoms. 
At second order in perturbation theory based on a double 
expansion of V so and a/R [14j, |15|, the SOI gives a spin- 
dependent shift °'^J ! g a of q T and directly shifts the en- 
ergy dispersion by — <7Ta2V so (a/R) cos 39, while the CIH 

only causes a valley-dependent q T shift T ^, F ™2 M ■ Here, 
the spin a = ±, the parameters, a\, ai, A relate to 
some unperturbed hopping amplitudes between carbon 
orbitals [16j , and 9 is the angle between the chiral vector 
and the zigzag direction along ai . Due to the hexagonal 
symmetry, this chiral angle is constrained to 0^0^30°, as 
calculated by fl=arctan 2 ^{f+n 2 f° r 0^n2^rii only. These 
corrections result in q T —> g CTT = r r ^ v + a ^ v s ° 1 f° r the low- 
est metallic n subband, e pTa — > s pGra —ahvF\fv 2 



<&t- 



ot A so2 , and V pTa -> V paTa = V a (K T + k)\ k ,., by 
setting A so i =a{V so a/R, A so2 =a 2 V so (a/R)cos39, and 
A cv —f3(a/ R) 2 cos 36*. Our generic Anderson Hamiltonian 
for a magnetic impurity coupled to the metallic CNT host 
then reads rl=% + H c , 

Rc ^ ^ P araC paTa C p(7Ta ^^V pa - T Q.Cp (TTa (i a -~{^li,C^j j , (1) 
p,a,T,a 

with the host DOS, p aT (e)=J2 p , a s ( e ~ Zpara), given by 



Paris) = p ( 
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where p =L/(hvp) and A CTT = |A CV + err A so i . Note that 
the BCS-like gap opens even in metallic CNTs. 

It is allowed to replace V paTa and £ paTa in the Hamil- 
tonian (1) by a proper constant coupling Vq and an effec- 
tive spectrum e poTOL , respectively, as long as the resultant 
effective DOS p(e) = J2 p T a S(e — s paTa ) satisfies p{e) = 
E p , T , a \VpaTa/V \ 2 S(s - e pr7Ta ) [17]. Remarkably, p(e) 



must be spin-independent due to the nonmagnetic nature 
of the CNT. By performing a Schrieffer- Wolff transforma- 
tion I18j , the system can then be readily mapped onto the 



Kondo model Rk— JS-s(O) which describes the antiferro- 
magnetic exchange interaction J=— 2VqU /[ed(sd+U)]>0 
of the impurity spin S with the host spin §(0) at the im- 
purity site 0, where Ed — Ed — Ep < is the impurity level 
measured from the Fermi energy Ep- 

The effective DOS p{e) seen by the impurity is essen- 
tially a renormalization of the bare CNT DOS, emerg- 
ing from the quantum interference between different hy- 
bridization paths the electrons can take to hop in and 
out of the impurity. Hereafter, we focus on a partic- 
ular class of impurity orbitals that hybridizes equally 
with the nearest carbon atoms on a given sublattice, i.e., 
V x j =V X - In this case, when the impurity is located on 
the S or C site, constructive interference renormalizes the 
CNT DOS as p sc (e)=E T ( £ + <JrA so 2 ) 2 PrTT (e) /{2Nt 2 ) by 
defining V = V b for S site or V = y/vj + V? for C site, 
whereas p T (e) — J2 T P<?t{z) / {'ZN) for T-sites adatoms 
where Vo=V a and the interference is absent. These DOS 
represent two distinct classes of hard-gap Kondo models 
promising for unconventional Kondo physics. 

(i) A half-filled (U = — 2s ( i) impurity coupled to the 
neutral armchair (9 = 30°) CNT, where since A SQ 2 = 
A cv = 0, no CIH effect exists and the SOI opens a gap of 
width, 2 A so i, at the Fermi level that exactly crosses the 
Dirac point. The system then exhibits strict p-h sym- 
metry, which prohibits the appearance of even powers of 
the local level £d in its renormalization by successively 
integrating out hig h energy states with energy ±A in the 



band edges [11 



19| . The lowest contribution arising from 
two-loop vertex renormalizations, up to the leading order 
in a double expansion of Vb and Ed, gives rise to the RG 
beta function (3(sd) = 4p(A)V r 2 e ( i/A. Consequently, the 
corresponding flow of the Kondo coupling J = — AVq /Id 
reads (3(J) = -Ap(k)V£ J / A. Solving this RG equation 
yields J (A) = J(A )exp{ - AV$ [p(A)/A 2 ] dA} , where 
Ao is the initial band cutoff. It is evident that consecu- 
tive RG transformations increase the effective coupling. 
However, as the scaling of the gapped Kondo model char- 
acterized by p sc (e) or p T (e), terminates at the gap edge 
A = A so i, </(A) flows to a finite value rather than infinity, 
signaling the absence of the strong coupling (SC) Kondo 
phase. Therefore, the impurity ground state is always 
a local moment, consistent with previous numerical RG 
calculations on the rectangularly gapped band f20j |. 

(ii) An mfinite-/7 impurity in the neutral armchair 
CNT. Here the p-h symmetry of the whole system is vi- 
olated, while the CNT bath remains p-h symmetric. RG 
transformations show that already at one-loop order, ver- 
tex renormalization of £d and J occurs as the band width 
is reduced 11, 0]. The resultant RG equation for the 
dimensionless Kondo coupling J = — 2 ( o(A)V r 2 /e ( 2 scales 
according to /3(J) = [lnp(A)]'AJ— J 2 . We solve the beta 



3 



function as 



J(A) 



p(A)J(Ao) 



p(A ) + J(A )/ j ( v o [p(A)/A]dA 



(3) 



where the denominator being vanishing or nonvanishing 
during scaling determines the impurity ground state. 
For T-site adatoms, the DOS p T (e) = p 8(|e| — 



= £o 

N ■ 



As 



A BO i)|e|/Ve 2 - A s 2 ol is of BCS-type, with p 
the scaling proceeds, the denominator of Eq. (3) van- 
ishes at the critical band width A c = T£ + A 2 ol /4T£ 
when 2T£ > A so i, directing the RG flow towards the 
SC Kondo fixed point. Here the Kondo temperature 
= Ao exp[— 1/ J(Ao)] is defined as a scaling invari- 
ant [ll| of the normal metallic model [realized by setting 
p(e) = p T (A ) ~ p for A > A so i]. By contrast, for 
2T° < A so i, or equivalently, J(A ) < [In (2A /A so i)]-\ 
the adatom flows to the unscreened LM state since the 
coupling J(A) already vanishes as the scaling enters 
into the gap region before it reaches the SC limit. We 
thus find a quantum-critical point of the impurity level 
£dc = —ZVqPo m 3^7 across which, upon lowering Ed, 
the impurity undergoes a QPT from a screened to an 
unscreened moment. The explicit R dependence of this 
phase boundary can be written as Idc = £ c i — 2VqP m 7T > 
with £cl = -2F 2 p o ln^-. 

The interference-induced additional scaling (e/t) 2 , im- 
posed on the effective host DOS p sc (e) — (e/t) 2 p T (e) 
for substitutional dopants or C-site adatoms, dramati- 
cally changes the RG flow of Eq. (3). It features a dif- 
ferent phase boundary at Idc — £ C 2 — /(-f ) separating 

the Kondo and LM phases, where e c2 = —VoPo~p~ an d 

/(£ ) = V Q % In by taking p s? (A ) ~ p § 

for A 3> A sol . For realistic parameters, this boundary is 
always much shallower than in the T-site case, reflecting 
a reduction of the Kondo regime by the interference. 

(iii) An infinite-?/ impurity coupled to the doped nan- 
otube with arbitrary chiralities. In this general case, 
the gaps in the two valley sectors of the CNT DOS 
are different, being centered at s = ±A S0 2 with width 
Wi = 2|A cv -A so i| and W 2 = 2(A cv +A so i), respectively. 
After summation over the two valleys, their overlap con- 
stitutes a net gap of width W — Wl + W2 — 2A so2 centered 
at e = min(A cv , A sol ) = e , in the effective DOS p T {e) 
and p sc (e), provided that W > 0. 

The presence of A so2 and A cv for the nonarmchair chi- 
rality, and the deviation of the Fermi level from the Dirac 
point in the doped nanotube, definitely break the p-h 
symmetry of the CNT baths. This renders the previous 
RG arguments insufficient because all vertex functions 
will develop structures on a scale Ef- The SB technique 
[2l| accounts for this complication by introducing an aux- 
iliary boson field to describe the empty state, together 
with a Lagrange multiplier A to exclude double occu- 
pancy in the impurity. The boson field is further con- 
densed to its saddle-point value r which minimizes the 



free energy, obeying ^ deIm[E(e)G(e)] = A, where 
= ilXde'r(e')(e-e' + i 0+)-\ T(e) = ^ 2 p(e), 
and the impurity Green's function G(e) = [e — Ed — A — 
r 2 E(e)] _1 . After A and r are self-consistently determined 
plj . the localized level Sd is renormalized to a Kondo 
resonance at Ed + A. This mean-field treatment correctly 
describes the Kondo fixed point for impurity levels deep 
below Ep y where charge fluctuations are frozen out. 
At the critical point, the Kondo temperature, defined 



as T K = v /(£ d + A) 2 +r 4 r 2 (e £i + A) [2l|, must vanish 



The saddle-point equation then yields the critical value 
of the impurity level, 



1 



~'dc 



1 



7T 7-A 6- E F 7T J E 



(4) 



This result is applicable to the doped case as well as to 
the nonarmchair chirality. Below we present numerical 
results for realistic CNT parameters [3,[l3|: ai=0.055, 
a 2 =0.217, /3=93.75meV, y so =6meV, and t=2.5eV. 

Figure 1 presents impurity phase diagrams in the (ed, 
R) plane, when the Fermi energy is tuned to the gap 
center, Ep — eq. The QPTs discussed here exist only if 
the CNT bath is gapped, such that p(Ep) vanishes ex- 
actly. Solving the inequation W > 0, we find a critical 
chiral angle 9o = | arccos(ai/a2) ~ 25.1° and an upper 
limit of CNT radius R = l3a/{V so a 2 ) w 17.7nm. For 
60 < ^ 30°, CNTs with arbitrary radius are always 
gapped (W > 0), resulting in that the impurity exhibits 
Kondo and LM ground states separated by transitions 
at Idc in the whole range of R [see, e.g., Figs. 1(a) and 
1(b) for the armchair case]. On the other hand, when 
^ 8 ^ 6q, the Kondo-LM transition can occur only 
for R < R a . Beyond this upper limit R ^ R , one has 
W ^ 0, leaving always a screened impurity state [see, 
e.g., Figs. 1(c) and 1(d) for the zigzag case]. The specific 
i?-dependence of the phase boundary is also very sen- 
sitive to explicit impurity positions, despite that p T (s) 
and P sc {e) feature the same gap structure. When the 
CNT radius increases, the Kondo regime of impurities 
on T (S or C) sites gradually widens (narrows), with 
the boundary eventually decreasing to — 00 (increasing 
to e c2 ) as R — > 00 for 9 > 9 chiralities [see, e.g., 
Figs. 1(a) and 1(b)] and R ->• R for 6 ^ 9 chiral- 
ities [see, e.g., Figs. 1(c) and 1(d)]. This sensitiveness 
stems from the quantum interference effect, which dra- 
matically changes the scaling behavior of p sc (e) out- 
side the gap region, as compared with p T (e). By not- 
ing e c i ~ 400e C 2 — — 10 8 /(1) < for the parameters 
used here, an interference-induced overall shrinking of 
the Kondo regime in the S /C configurations is also evi- 
dent, as already emphasized before. 

Gating the CNT host to tune its Fermi energy away 
from the gap center renders Ep closer to electronic states 
near one of the gap edges, in favor of screening the im- 
purity. Therefore, the SC fixed point can be reached for 
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FIG. 1: (Color online) Phase diagrams of Kondo-LM tran- 
sitions for impurities sitting on T [(a),(c)], S or C [(b), (d)] 
sites in the (n,n) armchair [(a),(b)] and (3n,0) zigzag [(c), (d)] 
CNTs, as a function of the impurity level Ed and discrete 
CNT radius R/a — y/Sn/2-K (armchair), 3n/2iv (zigzag), with 
b£Z. The high-energy cutoff Ao=0.5eV. 



smaller Kondo couplings (deeper impurity levels) , widen- 
ing the Kondo regime. This is confirmed by our calcu- 
lations shown in Fig. 2(a) for T-site adatoms [3], where 
the Kondo and LM phases are bounded by an arched 
borderline peaked at Ep = Eq. As Ep moves further 
out of the gap region, arbitrary small J > can al- 
ways drive the impurity into the Kondo phase due to 
p(Ep) 7^ 0. Interestingly, while the armchair CNT fea- 
tures a p-h symmetric phase diagram when it is tuned 
from hole doping (Ep < 0) to electron doping (Ep > 0), 
the arched phase boundary of nonarmchair CNTs always 
deviates from the hole-doped side because of e > 0, and 
can even fully enter into the electron-doped region for 
large CNT radius. The minimal radius R\ needed for ac- 
cessing this maximal p-h asymmetry can be determined 



a-2 cos 30 
Q1+Q2 cos 39 

14.1nm for zigzag CNTs. 



when 



by solving e > -^j- to obtain Ri = R 
9 < 6> , thereby Ri 

Obviously, for deep impurity levels, two consecutive 
QPTs occur whenever the Fermi energy sweeps over the 
two gap edges. These can be experimentally observed 
by scanning tunneling probes 0] which directly measure 
the impurity spectral densities, A(e) = — ^lmG(e). By 
placing the Fermi energy far away from the gap region 
[Fig. 2(b)], the smooth host DOS around Ep gives rise to 
conventional Kondo resonances in A(e). While the CNT 
chiralities are indistinguishable in these resonances, the 
interference inbuilt in the S /C configurations greatly nar- 
rows the resonances as compared with T sites, signaling 
a suppression of the Kondo effect which is not favorable 
for experimental observations, regardless of the existing 
universal scaling with the Kondo temperature [inset of 
Fig. 2(b)]. This scaling is violated when Ep is tuned to 
access the quantum critical region around the gap edges 
[Figs. 2(c) and 2(d)]. Specifically, although the scaling 
of different impurity positions persists to some extent, 
the zigzag CNT hosts two-peak Kondo resonances dis- 
tinct from the armchair one. We attribute the two-peak 
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FIG. 2: (Color online) (a) Phase diagrams for T-site adatoms 
upon variation of the Fermi energy Ef, with Ao = 0.5eV. (b)- 
(d) Spectral densities of the impurity with a deep level Sd = 
-25V t) 2 Po iri the (20,20) armchair and (36,0) zigzag CNTs 
for E f = eq- 0.88t (b), e - 0.55W (c), and e + 0.55W (d), 
with A = t. In (c) and (d), we set Vq = 0.1 W /[Txr 2 p(E F )] 
such that Tk ~ 0.1 W, and the curves corresponding to the 
armchair CNT are offset for clarity. 



structure to a distortion by the distinct DOS of nonarm- 
chair CNTs, which possess two singularities around each 
gap edge, arising from the two valley sectors of the bare 
DOS of Eq. (2). Finally, as Ep shifts into the gap region, 
the SB equations break down and the Kondo resonances 
immediately collapse into the featureless LM spectra. 

In conclusion, we have addressed the Kondo problem 
of magnetic impurities in CNTs, demonstrating the ex- 
istence of distinct QPTs in the impurity's ground state, 
which crucially depend on the characteristics of CNT and 
explicit impurity positions. 
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1) The parameters «i, ct2, and /3 

In the low-energy theory for carbon nanotubes, the ef- 
fects of spin-orbit interaction and surface curvature on 
7r electrons are well described in second-order perturba- 
tion theory [14,15]. The effects are equivalent to shift the 
dispersion relation by — <TTa 2 V so (a/ R) cos 38, to shift the 

perpendicular wave vector by aa\ ^° a R + t(3 a h ^° S j^ e , and 

to shift the parallel wave vector by t/3' a h ^^f ■ Assuming 
sufficiently long nanotubes, the last shift is irrelevant, we 
thus drop it. The explicit forms of the parameters, a±, 
a.2, and /3, appearing in the remaining terms are [14]: 



ai 



On 



V3e a (V£ + V£) 



18V s l 



V3v; p 



3(V n - V a Y 

uyv pp v pp j 

V (V v + V a ) 

pp \ pp 1 pp } 
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Here e s is the energy of the carbon s orbital relative 
to the p orbital energy. The latter (i.e., the on-site en- 
ergy of 7r electrons) is set to zero in our manuscript. V sp 
represents the unperturbed hopping amplitude between 
nearest-neighbor s and p orbitals. Vpp^ is the unper- 
turbed hopping amplitude between nearest-neighbor p 
orbitals, giving rise to the 7r(er) band. In this work, 
we use the parameter set: e s — — 8.9eV, V sp = 5.6eV, 
= -3.0eV, and = 5.0eV [R. Saito, G. Dressel- 
haus, and M. S. Dresselhaus, Physical Properties of Car- 
bon Nanotubes (Imperial College Press, London, 1998)], 
as it is also used in Ref. [14]. This give us, a\ ~ 0.055, 
c*2 — 0.217, and j3 ~ 93.75meV, to carry out numer- 
ical calculations. Using other sets of parameters [e.g., 
D. Tomanek and M.A. Schluter, Phys. Rev. Lett. 67, 
2331 (1991); J.W. Mintmire and C.T. White, Carbon 
33, 893 (1995)] does not change our numerical results 
qualitatively. 

2) Phase diagrams for S- and C-site impu- 
rities in the (e^, Ep) plane 



As shown in Fig. SI, the Kondo and LM phases of 
S/C impurities are also bounded by an arched borderline, 
showing features qualitatively same with T-site adatoms 
[see, Fig. 2(a) in the paper]. For example, the boundary 
is p-h symmetric for armchair nanotubes, but becomes p- 
h asymmetric for nonarmchair nanotubes. The minimal 
radius R\ derived in the paper for accessing the maximal 
p-h asymmetry also applies to this case. This is because 
p sc (s) and p T (e) share the same gap structure. They 
scale differently only outside the gap region due to the 
quantum interference effect. 

The effect of quantum interference is mainly reflected 
i) in the independence of the boundary (see Fig. 1 in 
the paper), ii) in the fact that the arched LM region 
of S/C impurities are much sharper than T adatoms 
[compare Fig. 2(a) in the paper with Fig. SI here], and 
iii) in the fact that for realistic nanotube parameters, the 
Kondo boundary of S/C impurities, Edc, is always much 
shallower than the boundary of T adatoms, signaling 
the reduction of Kondo regime by interference. 
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FIG. SI: Phase diagrams for substitutional dopants or C-site 
adatoms in armchair and zigzag nanotubes in the (ed, Ef) 
plane. The parameters used are the same as in Fig. 2(a). 



